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We investigate a potential of the higher multipole power spectra of the galaxy distribution in 
redshift space as a cosmological probe on halo scales. Based on the fact that a halo model explains 
well the multipole power spectra of the luminous red galaxy (LRG) sample in the Sloan Digital 
Sky Survey (SDSS), we focus our investigation on the random motions of the satellite LRGs that 
determine the higher multipole spectra at large wavenumbers. We show that our theoretical model 
fits the higher multipole spectra at large wavenumbers from N-body numerical simulations and we 
apply these results for testing the gravity theory and the velocity structure of galaxies on the halo 
scales. In this analysis, we use the multipole spectra Pa(Jz) and P&(k) on the small scales of the 
range of wavenumber 0.3 < fc/[/iMpc _1 ] < 0.6, which is in contrast to the usual method of testing 
gravity by targeting the linear growth rate on very large scales. We demonstrate that our method 
could be useful for testing gravity on the halo scales. 

PACS numbers: 


Redshift survey of galaxies is a promising way for exploring the nature the dark energy and testing 
gravity on the cosmological scales. Recent results of the baryon oscillation spectroscopic survey (BOSS) 
date release (DR) 11 of the Sloan Digital Sky Survey (SDSS) III have demonstrated the usefulness of 
redshift surveys [1, 2]. A possible tension in the cosmological parameters between the results by the 
Planck satellite and the BOSS is reported [3-6], which attracts interests of researchers. An interesting 
question that arises is whether this tension could be resolved in models where the gravity gets modified 
from its usual general relativistic form. 

The redshift-space distortion plays an important role in testing gravity [7, 8], which reflects the infor¬ 
mation of velocity of galaxies. One of the target of redshift surveys is a measurement of the redshift-space 
distortions in the linear regime of the density perturbations [9], which provides us with a chance of test¬ 
ing gravity through the linear growth rate. On the other hand, the Finger-of-God (FoG) effect is the 
redshift-space distortion in the nonlinear regime of density perturbations reflecting the random motion of 
galaxies. The primary purpose of the present paper is to investigate an effective method to evaluate the 
random velocity of galaxies in halos, which might provide us with a unique chance of testing gravity on 
halo scales. This can be achieved by precisely modeling the FoG effect on the basis of the halo model. 

In order to quantify the redshift-space distortions, the multipole power spectrum is useful, which is 
defined as a multipole coefficient of the multipole expansion of the anisotropic power spectrum (e.g., 
[6, 8, 10]). Recently, the authors of Ref. [11] found that a halo model describes well small-scale behavior 
of the higher multipole power spectra of the luminous red galaxy (LRG) sample of SDSS DR7. Based on 
this new finding, we consider a potential of measuring the velocity of satellite galaxies in halos and testing 
the gravity theory on the halo scales with the nrultipole power spectrum. The key of the method is the 
random motion of the satellite galaxies and their 1-dimensional velocity dispersion in a halo with mass 
M, for which we adopt a simple formula, 
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where (3 is a constant parameter, G is the Newton’s universal gravitational constant, and r v ; r is the virial 
radius defined by r v j r = (3M/47rp m (z)A v i r (z)) 1//3 , where p m (z) is the mean matter density and A v i r (z) 
is the density contrast of a halo, respectively, at the redshift z. We adopt A V i r = 265 at z = 0.3, for 
the sample corresponding to our LRG mock samples. We carefully check this velocity dispersion relation 
using the numerical simulations, as well as the validity of the theoretical model for the higher multipole 
power spectra. This theoretical model is compared with the SDSS LRG sample, and we put a useful 
constraint on the velocity dispersion and the gravitational constant on the halo scales. 

We here briefly review the multipole spectrum in a halo model according to Ref. [9, 11]. Following the 
general prescription of the halo approach [12-14], we write the anisotropic power spectrum in the redshift- 
space consisting of the 1-halo and 2-halo terms, PhRG(k, p) = P lh (k , p) + P 2h (k, p). We here consider the 
model which consists of the central galaxies and the satellite galaxies. We adopt the following expression 
(2) for the one halo term, 
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where we adopt the halo mass function dn/dM given by [15], n is the mean number density of LRGs 
given by n = f dM(dn/dM)NuoD(Ad), Nuou(AI) = (N cen ) + ( N sat ) is the halo occupation distribution 
(HOD), for which we adopt the following form [16], 
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with the error function erf(x), and cr 2 (M) and er 2 (M) are the velocity dispersion of the central LRGs and 
the satellite LRGs, respectively. Table 1 lists the HOD parameters matching the SDSS DR7 LRG catalog 
in Ref. [17]. We assume that the distribution of the satellite galaxies follows the NFW profile [18] and 
«NFw(fc) denotes the Fourier transform of truncated NFW profile [19]. Results of Ref. [20] support this 
assumption. We may assume that central LRGs reside near the halo center, thus their velocity difference 
relative to the host halo should be small (cf. [21]). On the other hand, satellite LRGs are off-centered and 
their random velocity should be the main source of the FoG effect. Here we assume 


cr 2 (M) = a 2 a 2 (Ad), 
a 2 (M)=a 2 s a 2 v (M), 
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where a c and a s are the constant parameters. 

In the previous paper [11], the 2-halo term was modeled with an analytic fitting formula from N-body 
simulations. However, in the present paper, we adopt a very simple treatment for the two halo term of 
the higher multipole power spectrum, because it is not trivial to construct a precise analytic model in 
redshift space which is applicable even at large wavenumbers. Using the mock catalogs corresponding the 
LRG sample is an alternative way to incorporate precise theoretical predictions for the two halo term. 
For this modeling, we adopt the results in the previous paper [9], which has constructed mock catalogs, 
corresponding to the SDSS LRG sample, and has investigated the behavior of the multipolc spectra. In 
the present paper, we use the following modeling for P 2h (k) and P§ h (k). The results in [9] demonstrate 
that the contribution from the two halo term to Pe(k ) is negligible, i.e., P|^(/c) ~ 0, and P 2h (k) is simply 
expressed as kP 2h (k) ~ lSf/iMpc^ 1 ] 2 , which we also adopt here. Because the contribution of the two halo 
term to P 2 (k) is rather large compared to that of Pi(k) and Pe(k) [11], it is not included in our analysis. 

We first demonstrate the validity of our theoretical model by comparing with the results of N-body sim¬ 
ulations. The simulations assume the spatially flat cold dark matter model with a cosmological constant, 
adopting D 0 = 0.273 and ag = 0.82. We run 10 realizations of N-body simulations using Gadget-2 code 
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TABLE I: HOD parameters of the LRG sample [17]. 


[22] with Gaussian initial condition. Each simulation has a side length of 600h -1 Mpc and the particle 
number of 800 3 (each particle mass is 2.8 x 10 1o /i _1 Mq). We use z = 0.3 snapshots and identify halos by 
Friends-of-Friends algorithm with a linking length of 0.2. Mock catalogs are constructed so that the bias 
and the HOD match the SDSS DR7 LRG catalog in Ref. [17]. The position of a central LRG is given by 
the potential minimum of the host halo and the velocity is given as the averaged velocity of all particles 
within the halos. We substitute randomly-picked up dark matter particles for satellite LRGs. In this 
analysis we constructed mock samples both with and without including the fiber collision effect [23]. We 
first make uncollided samples by removing one of the adjacent subhalos within 55 arcsec at z = 0.3 and 
randomly return a part of removed subhalos at 10 percent probability for the overlapped area of tiling 
where both spectrum of collided pairs can be measured. In our simulation, the central LRGs locate near 
the halo center, and their velocity is negligible. We assume no velocity bias for satellites. Thus, our mock 
catalogs should be understood as a c = 0 and a s = 1. 

Using the mock catalogs, we show the validity of our expression (1) for the velocity dispersion of satellite 
galaxies. The velocity dispersion of satellite galaxies in a halo has not been well understood, though there 
are a few works that investigate the velocity dispersion of LRGs [24, 25] . Recently, Guo et al. have studied 
the velocity bias of galaxies in the SDSS III CMASS sample in the context of a halo model [26]. Their 
results have implications for our results, as will be discussed below. 

Figure 1 compares the velocity dispersion er 2 (M) of satellites as a function of the host halo’s mass 
M. Here the cross symbols show the results of the N-body numerical simulation, while the curve shows 
(GM/2r v i r ) 1 / 2 , i.e., Eq. (1) with /3 = 1. This suggests that Eq. (1) with f3 = 1 reproduces well the relation 
between the velocity dispersion of satellite and the halo mass of our N-body simulations. 

The effect of the fiber collision, which misses galaxies located closely to each other, could be crucial 
in the analysis of the redshift-space clustering on small scales [27]. The fiber collision dominantly occurs 
for pairs in the same halo. In the previous work [11], the effect of the fiber collision is included by a 
multiplying factor reducing the satellite fraction. In the present paper, we adopt a similar prescription, 
for simplicity. Instead of introducing the satellite fraction, we float the HOD parameter Mi , which changes 
the satellite fraction, as a fitting parameter in our MCMC analysis. 

We compare the results of the numerical simulations and the theoretical multipole power spectra, with 
floating the two parameters /? and Mi. Especially, we here fixed a c = 0 and a s = 1 taking the consistency 
with our numerical simulations. Note that the HOD parameters other than Mi are fixed. In the MCMC 
analysis we only use P^k) and Pe(k) in the range of wavenumbers 0.3 < fc/f/iMpc -1 ] < 0.6 in order 
to reduce the influences from the uncertain contribution of the two halo term. Table II summarizes our 
results, where the best fitting values with one sigma statistical errors are presented for (A) simulation 
with the fiber collision (F.C.), (B) simulation without the fiber collision, and (C) LRG sample, (D) LRG 
sample with the two halo term modeled with the brightest LRG (BLRG) sample [11], from the left to the 
right column, respectively. The chi-squared and the degrees of freedom are also shown. In this table, the 
values within parentheses are the results for the data in the range of 0.2 < fc/f/iMpc^ 1 ] < 0.6. The left 
hand two columns of Figure 2 show the best fit curves of the HOD and the multipole power spectra for 
the simulations (A) and (B). These results demonstrate that our theoretical model reproduces those of 
the numerical simulations. 

We next apply our method to the multipole power spectra measured with the SDSS DR7 [28] . The DR7 
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FIG. 1: 1-dimensional velocity dispersion a v (M) as a function of halo mass. The crosses are from N-body 
simulation, while the curve is Eq. (1) with /3 = 1. 
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(B) Mock without F.C. 
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TABLE II: Results of our MCMC analysis with floating the 2 parameters /3 and Mi, where we fixed a c = 0 and 
a s = 1 and the other cosmological parameters. The best fitting values with one sigma statistical errors, the satellite 
fraction f dM(dn/dM)(N sa .t)/n, and the chi-squared along with the number of d.o.f are presented when fitted with 
(A) the simulation with the fiber collision (F.C.), (B) the simulation without the fiber collision, (C) LRG sample 
and (D) LRG sample with the two halo term modeled using the BLRG sample, from the left to the right column, 
respectively. The results are obtained using the data in the range of the wavenumbers 0.3 < fc/[/iMpc -1 ] < 0.6, and 
the values in the parenthesis are same but using the data in the range of the wavenumbers 0.2 < fc/[/iMpc _1 ] < 0.6. 


LRG sample is selected to cover a redshift range, 0.16 < z < 0.36, only in the northern cap in order to 
reduce systematic uncertainties and to match the analysis in Ref. [17]. Thus, the sky coverage is limited 
to 7189 deg. 2 and the total number of LRGs is 61 899. We adopt the same method for the measurement 
as that in Refs. [11, 29, 30], but with the fiducial cosmological background for the distant-redshift relation 
of the spatially flat ACDM cosmology with = 0.3. The right hand two columns in Table II list the 
results of the MCMC analysis with the LRG multipole spectra, whose difference comes from the modeling 
for the two halo term. The right two columns of Figure 2 show the best fitting curve and the data. 

The results of MCMC analysis with the SDSS LRG sample can be used for testing the gravitational 
constant on the halo scales. This is because the velocity dispersion in a modified gravity models could be 
written as = G e gM/2r v i r , where G e g is an effective gravitational constant. Regarding G e g = /3G, 

we may put a constraint on the effective gravitational constant from the SDSS LRG sample on the halo 
scales, /? = 1.70 ^q'® 5 from the column (C) in Table II, in which we adopted the same modeling for the two 
halo term as that of the mock catalogs. This value is rather larger than the prediction of the numerical 
simulations, although the error is not small. 

Though the contribution of the two halo term to P4 and Pq is rather small compared with the one halo 
term, but it might be influential to our results. As a check of our results, we model the contribution of 
the two halo term using the BLRG sample [11], Because the BLRG catalog roughly corresponds to the 
central galaxies catalog, then we may model the two halo term by computing the multipole spectrum of 
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the BLRG catalog. The column (D) of Table II represents the results f} = 1.35+1®. Compared with the 
case of the modeling for the two halo term from the numerical simulation, the value of /3 becomes small 
and /? = 1 is in the one-sigma error of the results. 

Let us discuss the reason why higher value of /3 is obtained from the analysis of the LRG sample. It could 
be a smoking gun of a modified gravity. For example, an f(R) gravity model has an effective gravitation 
constant G e s = 4G/3, as long as the chameleon mechanism does not work. However, we should discuss 
possible systematics that may lead to a larger value of /?. Because the satellite fraction is small, being 
around 6% of the total LRGs, the first term dominates the right hand side of Eq. (2). Then, taking the 
degeneracy in the central and the satellite galaxy velocities, we should understand that the constraint is 


+ + + 
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in the case (D) when we use the BLRG sample for modeling the two halo term. The results might be 
explained by a larger velocity dispersion of the central galaxy in the multiple system. Recently, Guo et 
al. have reported the velocity bias of galaxies in the SDSS CMASS samples [26]. The sample is different 
from ours, but they report that a c ~ 0.3. However, this value a c ~ 0.3 is rather small to explain our 
results, and a c ~ 0.6 is required within the general relativity /3 = 1. Other possible systematics is the 
modeling of the two halo term in Pe(k), as we obtained somewhat different values between (C) and (D) 
in Table II. More sophisticated simulations based on subhalo catalogs from N-body simulations could be 
necessary and useful. We checked that our treatment of the fiber collision effect works within the error. 
However, this effect could be more complicated and a more careful modeling of the fiber collision might 
be necessary. 

In summary, we have investigated the potential of the higher multipole power spectra of the galaxy 
distribution in redshift-space. This method is based on the recent finding that a halo model accounts 
well for the behavior of the multipole power spectrum of LRGs on small scales. Our method uses the 
data of the spectrum on small scales 0.3 < fc/[/iMpc+ < 0.6. This is quite in contrast to the usual 
method of testing gravity by measuring the linear growth rate on very large scales. Our method is based 
on the fact that one halo term makes a dominant contribution to the higher multipole power spectra at 
large wavenumbers, which reflects the random motions of the satellite galaxies. We carefully investigated 
the relation between the velocity dispersion of the random motions of satellite galaxies and the host halo 
mass on the basis of the mock catalogs from N-body simulations. The validity of our theoretical model 
for the higher multipole power spectrum is tested using the results of the mock catalogs. By confronting 
our theoretical model and the observed multipole spectra of the SDSS LRG samples, we obtained a value 
for an effective gravitational constant somewhat larger than that predicted by the numerical simulations. 
This could be a smoking gun of the modified gravity. However, we might need to check our theoretical 
model for the two halo term and the fiber collision effect more carefully. 

K. Y. thanks the workshop, APC-YITP collaboration: Mini-Workshop on Gravitation and Cosmology, 
which was held at YITP Kyoto universe for a useful chance for discussions on the topic of the present 
paper. This work is supported by a research support program of Hiroshima University. The research by 
C.H. is supported in part by Grant-in-Aid for Scientific researcher of Japanese Ministry of Education, 
Culture, Sports, Science, and Technology (No. 24740160). 


[1] L. Anderson, et al., Mon. Not. Roy. Astron. Soc. 441, 24 (2014) 

[2] W. J. Percival, et al., Mon. Not. Roy. Astron. Soc. 439, 2531 (2014) 

[3] L. Samushia, et al., Mon. Not. Roy. Astron. Soc. 439, 3504 (2014) 

[4] C-H. Chuang, et al., arXiv: 1312.4889 

[5] A. G. Sanchez, et al., Mon. Not. Roy. Astron. Soc. 440, 2692 (2014) 

[6] F. Beutler, et al., Mon. Not. Roy. Astron. Soc. 443, 1065 (2014) 

[7] L. Guzzo, et al., Nature 451, 541 (2008). 

[8] K. Yamamoto, T. Sato, G. Huetsi, Prog. Theor. Phys. 120, 609 (2008) 




6 


[9] C. Hikage, Mon. Not. Roy. Astron. Soc. 441, L21 (2014) 

[10] K. Yamamoto, et al., Publ. Astron. Soc. Japan 58, 93 (2006) 

[11] C. Hikage, and K. Yamamoto, J. Cosmol. Astropart. Phys. 1308, 019 (2013) 

[12] U. Seljak, Mon. Not. Roy. Astron. Soc. 318, 203 (2000) 

[13] M. White, Mon. Not. Roy. Astron. Soc. 321, 1 (2001) 

[14] A. Cooray and R. Sheth, Phys. Rep. 372, 1 (2002) 

[15] R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 308, 119 (1999) 

[16] Z. Zheng et al., Astrophys. J. 633, 791 (2005) 

[17] B. A. Reid and D. N. Spergel, Astrophys. J. 698, 143 (2009) 

[18] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 490, 493 (1997) 

[19] R. Scoccimarro, R. K. Sheth, L. Hui, and B. Jain, Astrophys. J. 546, 20 (2001) 

[20] H. Guo, et al., Mon. Not. Roy. Astron. Soc. 441, 2398 (2014) 

[21] C. Hikage, et al. Mon. Not. Roy. Astron. Soc. 435, 2345 (2005) 

[22] V. Springel, Mon. Not. Roy. Astron. Soc. 364, 1105 (2005) 

[23] M. R. Blanton et al., Astron. J. 125, 2276 (2003) 

[24] C. Hikage, M. Takada, D. N. Spergel, Mon. Not. Roy. Astron. Soc. 419, 3457 (2012) 

[25] S. Masaki, et al., Mon. Not. Roy. Astron. Soc. 433, 3506 (2013) 

[26] H. Guo, et al., Mon. Not. Roy. Astron. Soc. 446, 578 (2015) 

[27] H. Guo, I. Zehav, Z. Zheng, Astrophys. J. 756, 127 (2012) 

[28] K. N. Abazajian, et al., Astrophys. J. Suppl. 182, 543 (2009) 

[29] A. Oka, et al., Mon. Not. Roy. Astron. Soc. 439, 2515 (2014) 

[30] K. Yamamoto, et al., Phys. Rev. D 81, 103517 (2010) 


7 


(A) Mock with F.C. (B) Mock without F.C. 
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FIG. 2: Top panels show the HOD N S (M) (dotted lines) and N C (M) (solid lines), middle panels kP^{k), and 
bottom panels kPe(k). In each panel, the curves correspond to the best fitting models. Each column shows the 
results of (A) mock with fiber collision, (B) mock without fiber collision, (C) LRG sample, and (D) LRG sample 
subtracted the two halo term modeled using the BLRG sample, from left to right, respectively. 



















